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The characteristics of a rotating stall of an impeller and diffuser and the evolution of a vortex generated at the 
diffuser leading-edge (i.e., the leading-edge vortex (LEV)) in a centrifugal compressor were investigated by expe- 
riments and numerical analysis. The results of the experiments revealed that both the impeller and diffuser rotat- 
ing stalls occurred at 55 and 25 Hz during off-design flow operation. For both, stall cells existed only on the 
shroud side of the flow passages, which is very close to the source location of the LEV. According to the CFD re- 
sults, the LEV is made up of multiple vortices. The LEV is a combination of a separated vortex near the lead- 
ing-edge and a longitudinal vortex generated by the extended tip-leakage flow from the impeller. Therefore, the 
LEV is generated by the accumulation of vorticity caused by the velocity gradient of the impeller discharge flow. 
In partial-flow operation, the spanwise extent and the position of the LEV origin are temporarily transmuted. The 
LEV develops with a drop in the velocity in the diffuser passage and forms a significant blockage within the dif- 


fuser passage. Therefore, the LEV may be regarded as being one of the causes of a diffuser stall in a centrifugal 


compressor. 
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Introduction 


Centrifugal compressors equipped with vaned diffu- 
sers are frequently employed in industry due to their high 
pressure-rise characteristics. However, they have some 
shortcomings, such as a narrower steady operating range 
than those equipped with vaneless diffusers. Thus, un- 
steady phenomena such as surges and rotating stalls are 
likely to occur at off-design points. Turbochargers are 
often used under off-design conditions, but centrifugal 
compressors should not be subjected to such operating 
conditions. Therefore, clarifying the unsteady phenomena 
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that could cause a serious accident is important to im- 
proving not only the performance of the compressor but 
also the safety at the off-design points. 

Many studies have examined the features of the rotat- 
ing stalls that occur in centrifugal compressors with 
vaneless diffusers. Senoo et al. (1977, 1978) were the 
first to investigate the relationship between the reverse 
flow and the inception of a rotating stall in the vaneless 
diffuser of a centrifugal blower, and they proposed crite- 
ria for the inception of a rotating stall in a vaneless dif- 
fuser [1,2]. Haupt et al. (1988) investigated the pressure 
fluctuation that occurred between the impeller inlet and 
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Nomenclature 
B diffuser passage height (m) 
D diameter (m) 
Le chord of diffuser vane (m) 
La distance between impeller exit and diffuser 
leading-edge (m) 
N rotational speed (min’') 
Q volume flow rate(m*/s) 
P static pressure (Pa) 
P, total pressure (Pa) 
Ur circumferential velocity at impeller exit (m/s) 
Vr radial velocity (m/s) 
Greek Letters 
$ flow coefficient 
Wr total pressure-rise coefficient 


the diffuser outlet, and found that it was caused by peri- 
odic changes in the flow conditions near the shroud side 
of the impeller inlet [3]. Iwakiri et al. (2009) reported 
that the tornado-type separation vortex caused by 
full-blade leading-edge separation dominated the flow 
field of a developed rotating stall condition. There have 
been a few reports about the spike and modal stall incep- 
tions in a centrifugal compressor with a vaned diffuser 
[4]. Spakovszky et al. (2007) reported that the path to 
instability is governed by either the spike or modal de- 
pending on whether the slope of the pressure-rise charac- 
teristics in the semi-vaneless space becomes zero before 
the peak of the overall diffuser static pressure rise cha- 
racteristics is reached [5]. 

However, very few researchers have reported the un- 
steady phenomena in centrifugal compressors equipped 
with vaned diffusers. The causes of a rotating stall, espe- 
cially the relationship between the diffuser stall and an 
unsteady flow field, are not yet fully understood. The 
authors reported on the unsteady behavior and evolution 
of the leading-edge vortex (LEV) in off-design operation, 
which was proposed to be a possible cause of a diffuser 
stall in a centrifugal compressor [6, 7]. 

The present study focused on the impeller and diffuser 
rotating stalls that occur in the vaned diffuser of a centri- 
fugal compressor during off-design operation. Specifi- 
cally, we examined the relationship between the diffuser 
stall and unsteady behavior of the LEV. The flow fields 
in the diffuser passage were visualized by using the 
oil-film method. Both diffuser and impeller stalls were 
detected in the partial flow operation by the application 
of unsteady pressure and velocity measurements. In the 
computational fluid dynamics (CFD) analysis, detached 
eddy simulation (DES) was used to establish numerical 
schemes to investigate the vortex structure in detail. Ac- 
cording to the CFD analysis, the LEV was a combination 


J. Therm. Sci., Vol.25, No.1, 2016 


Subscripts 

1 impeller inlet 

2 impeller outlet 

3 diffuser leading-edge 

4 diffuser trailing-edge 
Abbreviations 

BPF blade-passing frequency 
D.E. diffuser exit 

D.I. diffuser inlet 

LI. impeller inlet 

LEV leading-edge vortex 
ODV original wedge-type diffuser vane 
VL vaneless diffuser 


of a separated vortex near the leading-edge and a longi- 
tudinal vortex generated by the extended tip-leakage flow 
from the impeller. Additionally, the LEV formed a sig- 
nificant flow blockage within the diffuser passage when 
the velocity decreased within the passage. The results 
showed that a diffuser stall may be caused by the evolu- 
tion of the LEV. 


Experimental Apparatus and Procedure 


Experimental Apparatus 


The tested compressor was a centrifugal compressor 
having a design that was based on a turbocharger for ma- 
rine diesel engines. The geometric shape and dimensions 
of the tested compressor are listed in Table 1, and a sec- 
tional view of the compressor is shown in Fig. 1. The 


Table 1 Dimensions of Tested Compressor 
Tested Centrifugal Compressor 


Rotational Speed N 6000 min 
Rate Flow Mass G 1.64 kg/s 
Pressure Ratio P;/Po 1.1 
Impeller 
Number of Blades Z 14 
(Main + Splitter) (7+7) 
Inlet Diameter D; 248 mm 
Outlet Diameter D, 328 mm 
Exit Blade Width By 26.14 mm 
Diffuser 
Blade Shape Wedge 
Number of Vanes V 15 
Leading-edge Diameter D; 360 mm 
Trailing-edge Diameter Ds 559 mm 
Diffuser Width By 26.14 mm 
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Revolution Measurement System 


Pulse Generator 


15 Diffuser Vanes 


Oil-Film Method 


E Pressure Transducer 


Hot-wire 
Anemometer 


Fig. 1 Experimental Apparatus and Measuring System. 


tested impeller was an open-shroud type with seven main 
and splitter blades, and a specific speed of 401.8. The 
outlet angle of the impeller was 56°. Its inlet and outlet 
diameters were 248 and 328 mm, respectively. The shroud 
clearance of the impeller was about 0.2% of the impeller 
inlet blade width, B; = 105.6 mm. 

Two types of diffusers were used in the experiments: 
vaneless and vaned. The vaned diffuser had 15 diffuser 
vanes installed in the diffuser passage. The distance be- 
tween the impeller exit and the diffuser leading-edge Lig 
was set to 16 mm. The vanes were located between two 
parallel diffuser walls with heights B, of 26.14 mm. The 
original wedge-type diffuser vane, abbreviated to ODV, 
was used in the experiments. The vanes were used in the 
experiment by attaching them to the shroud surface of the 
channel diffuser. 


Measuring Method 


The rotational speed N of the tested compressor was 
limited to 6000 min”. The operation point was set by a 
butterfly valve installed at the outlet duct. The mass flow 
rates were calculated by using an orifice flow meter and a 
thermo-couple installed at the outlet duct. The pressure in 
the compressor was measured using a differential pres- 


sure transducer mounted at the entrance to the outlet duct. 


The flow field in the diffuser passage was visualized by 
using the oil-film method. 

The pressure and velocity measuring systems are also 
shown in Fig. 1. The pressure fluctuation on the shroud 
surface of the compressor was measured by using pres- 
sure transducers (Kulite, XCQ-062-25A) mounted in 
three locations. The pressure transducers were mounted 
flush with the shroud surface of the compressor. The 
pressure was measured at the impeller inlet, impeller exit, 
and diffuser outlet. The unsteady velocity spanwise dis- 
tribution of the impeller-discharge flow was measured at 
D = 339 mm, which was located between the impeller 


exit and the diffuser leading-edge, with an I-type hot-wire 
anemometer (DANTEC, 55R01). The impeller-discharge 
flow traversed from the hub to the shroud side every 5% 
of the diffuser width. 


Computational Procedure 


Governing Equations 


The three-dimensional unsteady flow field was calcu- 
lated by means of numerical simulation to investigate the 
flow fields in the impeller and diffuser passage. The 
simulations were carried out using in-house CFD code by 
solving the governing equation of continuity, a 3-D com- 
pressible Reynolds-averaged Navier-Stokes equation, an 
energy equation, and an equation for the state of an ideal 
gas. The inviscid flux was evaluated using flux difference 
splitting (FDS) [8], and extended to a high order using 
the monotonic upwind scheme for conservative laws 
(MUSCL) interpolation [9]. The viscous flux was discre- 
tized as the second-order central difference with Gauss’ 
theorem. The matrix free Gauss-Seidel (MFGS) implicit 
algorithm was used for the time integration [10]. In the 
turbulence model, detached eddy simulation (DES) was 
adopted; this is a hybrid scheme involving LES and 
RANS [11]. In this study, DES was based on the SST k-a@ 
turbulence model [12]. This model includes a depend- 
ency on the local turbulence length. Both the Coriolis and 
centrifugal forces were regarded as being inertial force 
terms in the relative coordinate systems. 


Computational Domains 


An overview of computational grids used in the nu- 
merical simulation is illustrated in Fig. 2. A composite 
grid system with structured H-type grids was used to 
simulate the flow field in the impeller and diffuser pas- 
sages. The computational domain was divided into three 
regions: the moving impeller region, the stationary dif- 
fuser region, and the moving tip clearance region. 


Fig. 2. Overview of Computational Domain. 
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In this study, we first analyzed one blade-to-blade im- 
peller and three diffuser passages to investigate the de- 
tails of the mechanism whereby a LEV is generated. The 
impeller (including the clearance region) and diffuser 
system had 3,523,368 and 5,177,361 cells, respectively. 
The total grid system had 8,700,729 cells. Additionally, 
we also conducted a full annulus analysis to investigate 
the impeller and diffuser rotating stall. The computational 
grids incorporated 7 main and splitter blades and 15 dif- 
fuser vanes. In total, the computational grids consisted of 
51,072,178 cells. 


Boundary Conditions 


At the inflow boundary, the total pressure and total 
temperature were fixed using the Riemann invariant of 
the one-dimensional characteristic waves to establish a 
non-reflecting condition [13]. At the outflow boundary, in 
addition to the Riemann invariant, a throttle resistance 
with a load coefficient was assumed between the bound- 
ary and the external atmosphere. Therefore, the static 
pressure at the outlet boundary was related to the mass 
flow rate through the boundary. 

Across the sliding boundary separating the moving 
impeller and the stationary diffuser frames, the latest data 
on one side was interpolated by the data for the opposite 
side by using a sliding mesh for unsteady simulation. 
Non-slip and adiabatic conditions were adopted for the 
wall conditions. The chimera grid method was adopted to 
interpolate the conservation values between the moving 
tip clearance and impeller regions [14]. 


Results and Discussion 


Characteristics of Tested Compressor 


The measured results for the compressor performance 
are shown in Fig. 3. The figure also shows the operating 
points obtained by DES analysis for a partial annulus. 
The flow and total pressure-rise coefficients are defined 


as follows: 
4 m’ D?B, N’ á pr,D?N’/2 @ 

The ODV, indicated by the red circles in the figure, 
was able to achieve high pressure-rise characteristics 
relative to those of the vaneless diffuser (VL), indicated 
by the black squares. However, the steady operation range 
became narrower, and the noise level increased remarka- 
bly [15,16]. 

In the designed flow operation (¢ = 0.24), the results 
of unsteady CFD analysis using DES were in good 
agreement with the measured results obtained using the 
ODV. However, the simulation predicted a higher total 
pressure rise than was actually measured at an off-design 
point. One reason for this is thought to be that a numeri- 
cal analysis of a partial flow passage cannot predict the 
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unsteadiness of the flow field. 


Mechanism Whereby Leading-edge Vortex is Formed 


The visualization results for the flow field on the 
shroud and suction surfaces of the ODV, as obtained us- 
ing the oil-film method are shown in Fig. 4. The LEV 
was observed as being stationary on the suction side of 
the diffuser leading-edge, as indicated by the white cir- 
cles in the figure. The size of the LEV developed to more 
than 80% of the diffuser leading-edge height at off-design 
operation of ø = 0.14. According to the visualization re- 
sults obtained using the colored oil film method, the flow 
field at the shroud surface between the impeller exit and 
the diffuser leading-edge was reversed at flow coeffi- 
cients of both ø= 0.24 and ¢= 0.14. A branching stream- 
line was observed near the diffuser leading-edge during 
design operation of ¢= 0.24. However, it moved into the 
diffuser inner passage when the operation point shifted to 
¢=0.14. The reverse flow appeared to be much stronger, 
such that the reverse flow region became much larger. 
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Fig. 3 Compressor Performance 
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(ii) Shroud Wall (¢ = 0.14) 
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Fig. 4 Visualization of Diffuser Passage Flow and Leading- 
edge Vortex (LEV) (ODV, ®= 0.24 and 0.14). 
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Fig. 5 shows the visualization results for the LEV ob- 
tained both experimentally (N = 6000 min`’) and by nu- 
merical analysis (N = 22,000 min"). The design rotation- 
al speed of the tested compressor was 22,000 min’ for 
practical use. The LEV had the same form both at N = 
22,000 min! and N = 6000 min’. Therefore, the experi- 
ments and numerical analysis were carried out at a low 
rotational speed of 6000 min” to enable easier visualiza- 
tion and minimize disturbances caused by unsteady phe- 
nomena such as passage shock waves. The results of the 
low-speed experiment and visualization may enable us to 
elucidate the unsteady behavior of the LEV and can be 
generalized not only for high-pressure compressors but 
also for other low-pressure fluid machinery such as cen- 
trifugal fans and blowers. 


i spe LAY 
Shroud < 


(i) Experiment (N = 6000 min”) 


Gi) CFD (N = 22,000 min) 


Fig. 5 Visualization of Leading-edge Vortex (LEV) for Rota- 
tional Speed Change(ODV, ®= 0.24). 


According to the numerical visualization of the vortex 
by Q-definition [17], colored by non-dimensional helicity, 
the LEV was made up of two longitudinal vortices at the 
design point (Fig. 6). From previous results, it was found 
that the incidence angle of the diffuser vane was ex- 
tremely positive near the shroud side of the diffuser 
leading-edge, and a transverse vortex was generated near 
the shroud side of the diffuser suction surface [6]. Fur- 
thermore, the radial velocity on the shroud side was 
much lower than that on the hub side near the impeller 
exit, and the flow field on the shroud side near the dif- 
fuser leading-edge was reversed by the high adverse 
pressure gradient. Therefore, vorticity was continuously 
generated, and two longitudinal vortices A and B were 
produced near the shroud side of the diffuser vane. Lon- 
gitudinal vortex A was dragged by the transverse vortex 
and collided with the diffuser suction surface. On the 
other hand, longitudinal vortex B, which was not dragged 
to the diffuser suction surface, moved toward the pres- 
sure surface of the adjacent diffuser vane. 

Fig. 7 shows the instantaneous vortical structure in the 
impeller-fixed relative flow field within the impeller and 
diffuser passage. A tip-leakage vortex, identified by the 
blue trace, was produced near the leading-edge of the 
impeller blade and extended to the diffuser passage, as 
shown in Fig.7 (i). This vortex was also visible in Fig. 7 
(ii) as a red trace and forms a part of the LEV by im- 
pinging on the leading-edge of the diffuser vane. There- 


fore, two longitudinal vortices, shown in Fig. 6, were 
produced by the tip-leakage vortices from the impeller. 
From these results, the LEV was determined to be a 
combination of a separated vortex near the leading-edge 
and a longitudinal vortex generated by the extended 
tip-leakage flow from the impeller. 


Jl 


Transverse 


Vortex 


Fig. 6 Structure of Leading-edge Vortex (LEV) (ODV, ® = 
0.24). 
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Fig. 7 Origin of Leading-edge Vortex (LEV) (ODV, ®= 0.24). 


Additionally, the instantaneous vorticity distributions 
near the leading-edge of the diffuser vane at an off- 
design point are shown in Fig. 8. In the flow fields at 
time t = 0, the LEV, indicated by the yellow circles, ex- 
isted on the shroud side of the diffuser suction surface. A 
tip-leakage vortex from the impeller, as shown by the 
white circles, was also observed near the leading-edge of 
the diffuser vane. Then, the LEV was extended to the 
adjacent diffuser vane by the impeller discharge flow and 
developed in combination with the tip-leakage vortex at 
time t = 0.5t,,. A new tip-leakage vortex, indicated by the 
red circles, was also produced by the impeller and moved 
toward the leading-edge of the diffuser vane. Therefore, 
the LEV was enlarged by the accumulation of the vortic- 
ity caused by the tip-leakage vortex. 


Characteristics of Impeller and Diffuser Stall 


To investigate the unsteady characteristics of the im- 
peller and diffuser stall, the surface static pressure fluctu- 
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ations were simultaneously measured using high-response 
pressure transducers mounted flush with the shroud sur- 
face at three different locations (Fig. 1); the impeller inlet 
(I.1.), diffuser inlet (D.I.), and diffuser exit (D.E.). The 
simultaneous pressure measurements along the circum- 
ferential direction were conducted using pressure trans- 
ducers installed at D.I.1 and D.I.2. The angle between 
D.I.1 and D.I.2 was 24°. 


Vortex 
From Tip 


Vortex 
From Tip 


Diffuser 
(i) Diffuser Passages (ii) Diffuser Suction Surface 
Fig. 8 Enhanced Structure of Leading-edge Vortex (LEV) 
(ODV, ®= 0.14). 


The coherence functions of the pressure fluctuation 
obtained at the four positions (I.I., D.I.1, D.I.2, and D.E.) 
during partial flow operation using ODV are shown in 
Fig. 9. Two remarkable peaks can be seen at 25 and 55 
Hz in the coherence functions measured at (D.I.1, D.I.2) 
and (D.I., D.E.), excluding the peaks for the fundamental 
and half blade-passing frequency (BPF) components. 
However, in the coherence function measured at (D.L, 
LI), the 25-Hz peak was attenuated, such that only the 
55-Hz peak was observed. As shown in Fig. 10, the 55- 
Hz peak was also recognizable in the coherence function 
measured at (D.I., I.I) at the off-design point of ø= 0.13 
when the VL was installed. Therefore, the 25-Hz peak 
may occur only within the diffuser passages because of 
the local unsteady phenomena, and the 55-Hz peak may 
originate within the impeller passages. The frequency of 
the two peaks shown in Fig. 9 moved linearly when the 
compressor rotational speed changed [7]. 

From the above measured results, the pressure fluctua- 
tions measured at 25 and 55 Hz were found to be caused 
by rotating stall within the diffuser and impeller passages, 
respectively. The rotating speeds of the stall cells within 
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the diffuser and impeller passages were easily estimated 
to be 25% and 55% of the impeller rotational speed. 

The impeller discharge velocity was measured with a 
hot-wire anemometer to investigate the spanwise distri- 
bution of the velocity fluctuation and the stall cell beha- 
vior. The spatial distribution of the FFT results in the 
spanwise direction is shown in Fig. 11. The horizontal 
axis is the frequency, and the vertical axis is the spanwise 
direction measured from the shroud surface. The contour 
colors express the magnitude of the velocity fluctuation 
normalized by the averaged circumferential velocity 
measured at each span. 

The intensities of both the diffuser and impeller stall 
fluctuations showed large magnitudes at the location be- 
tween the shroud side and mid-span (Fig. 11). The inten- 
sity and scale of the diffuser stall fluctuation were much 
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Fig.9 Coherence of Pressure Fluctuation (ODV, = 0.14). 
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Fig. 10 Coherence of Pressure Fluctuation (VL, ®= 0.13). 
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larger than those of the impeller stall fluctuation. Because 
both the impeller and diffuser stalls occurred only on the 
shroud side in the spanwise direction, which is very close 
to the source location of the LEV, the unsteady behavior 
of the LEV may play an important role in the inception of 
the rotating stall. Therefore, the relationship between the 
LEV and the diffuser stall inception process was investi- 
gated. 


Velocity Fluctuation % 


100 


Span from Shroud H % 


0 50 10? 
Frequency Hz 


Fig. 11 Distribution of Velocity Fluctuation in Spanwise Di- 
rection (ODV, ®= 0.14). 


Unsteady Behavior of Leading-edge Vortex 


From the results shown in Fig. 4, the flow field on the 
shroud side can be seen to change direction toward the 
adjacent vane in the diffuser. Then, the scale of the LEV 
increases in the case of partial flow operation. Thus, the 
LEV may develop in the direction of the adjacent vane 
and block the main stream in the diffuser passage. This 
enlargement of the LEV was assumed to increase the 
backflow rate and reverse flow region on the shroud wall. 
The behavior of the LEV can be observed in the lower 
part of Fig. 4, where it blocks the flow on the shroud side 
in the diffuser passages. The evolution of the flow 
blockage was considered to be caused by the diffuser loss 
and blade stalling in the diffuser passage. 

According to the CFD results for the partial flow pas- 
sages, the velocity measured at the diffuser leading-edge 
is shown in Fig. 12. The velocity was calculated from the 
mass flow passing through the throat of the diffuser vane. 
LP (red line) and RP (blue line) in the figure indicate the 
left and right sides of the diffuser passages, respectively, 
as shown in Fig. 13. The passages were beside each other. 
In the designed flow operation (ø = 0.24), the velocities 
of LP and RP agreed with each other. The instantaneous 
flow fields at times tg), ta2, and taz in Fig. 12 are shown 
schematically in Fig. 13. In these figures, the LEV beha- 
vior is visualized with the Q-definition and colored by 
the non-dimensional helicity. In the flow fields at time t = 
ta, the LEV generated from diffuser A, especially longi- 
tudinal vortex B in Fig. 7, moved toward the impeller 
passages because the origin of the longitudinal vortex B 
was the tip-leakage vortex from the impeller. Then, the 


LEV moved toward the leading-edge of diffuser B as a 
result of the impeller discharge flow at time t = tg. For 
the suction side of diffuser A, the center of the LEV at 
time t = tg; moved downstream, and the center of a new 
LEV was generated. Then, the LEV became attached to 
the next diffuser leading-edge. Therefore, the center of 
the LEV at time t = ta; on the suction side of diffuser A 
disappeared. This implies that the tip-leakage vortex syn- 
chronized with the impeller rotation and impinged on the 
leading-edge of each diffuser vane. From the above re- 
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Fig. 12 Fluctuation of Diffuser Inlet Velocity (ODV) 
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Fig. 13 Unsteady Behavior of Leading-edge Vortex (ODV, 
p= 0.24). 
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sults, the propagation of the LEV in the diffuser passages 
influenced both the impeller and diffuser passages near 
the shroud side despite the design operation. However, 
the size of the LEV did not change significantly, and the 
LEV was comparatively stable. Therefore, the flow fields 
in the diffuser passages also appeared to be stable. 

On the other hand, a large difference was observed at 
the off-design point (Fig. 14). In particular, the velocity 
of LP was much lower than that of RP, and a diffuser stall 
occurred in LP. The size of the LEV at time t = t,; was 
greater than that at the designed flow operation, and the 
size of the LEV and position of its origin were tempora- 
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Fig. 14 


Visualization of Vortical Structure by Full Annulus 
Analysis 

Fig. 15 shows the instantaneous vortical structure for 
design and off-design operation using ODV as obtained 
by full annulus analysis. The vortical structure was visu- 
alized with the Q-definition and colored by the 
non-dimensional helicity. The LEV was generated on the 
suction side of the diffuser vanes at the design point as 
predicted by partial passage analysis, as shown in Fig. 15 
(1). The tip-leakage vortex was also recognized and con- 
sidered to influence both the impeller and diffuser pas- 
sages. Additionally, Fig. 15 (ii) shows that the LEV was 
enlarged and extended to the adjacent diffuser vane. The 
spillage of the LEV from one diffuser passage to another 
was considered to be related to the diffuser rotating stall 
that may occur in multi-passages. However, the mechan- 
ism of impeller and diffuser stalls was not fully unders- 
tood. Therefore, further investigation of impeller and 
diffuser stalls will be conducted by the application of full 
annulus analysis in our future work. 
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rily transmuted. The LEV was expanded over the entire 
span by the accumulation of the vorticity caused by the 
velocity gradient resulting from the impeller discharge 
flow at times t,- t,3. At time t,3- tp4, the size of the LEV 
fell as the velocity in the diffuser passage increased. The 
above results indicate that the LEV blocked the flow in 
the diffuser passages, especially the longitudinal vortex B 
shown in Fig. 7, when the mass flow rate was decreased, 
and this was assumed to be the cause of the diffuser stall 
in the centrifugal compressor diffuser. Therefore, con- 
trolling the LEV is important for improving the com- 
pressor performance during off-design operation [6]. 
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Fig. 15 Vortical Structure of full annulus analysis (ODV) 
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Conclusions 


The characteristics of impeller and diffuser stalls and 
diffuser leading-edge vortex in a centrifugal compressor 
were investigated by experiment and through CFD anal- 
ysis. 

The results can be summarized as follows: 

(1) The LEV was observed on the shroud side of the 
diffuser suction side near the leading-edge. The LEV 
could be observed even at the design point and devel- 
oped into a flow blockage in the diffuser passage during 
off-design flow operation. 

(2) The LEV consisted of two longitudinal vortices A 
and B. Longitudinal vortex A was dragged by the trans- 
verse vortex and collided with the diffuser suction sur- 
face. On the other hand, longitudinal vortex B was not 
dragged to the diffuser suction surface and instead 
moved toward the pressure surface of the adjacent dif- 
fuser vane. 

(3) A tip-leakage vortex was generated near the lead- 
ing-edge of the impeller blade and forms a part of the 
LEV by impinging on the leading-edge of the diffuser 
vane. Additionally, the LEV was enlarged by the accu- 
mulation of the vorticity caused by the tip-leakage vor- 
tex. 

(4) In the tested compressor, both the impeller and 
diffuser stalls occurred at 55 and 25 Hz during off-design 
flow operation. Both stall cells were present only on the 
shroud side of the flow passage, where they rotated. The 
diffuser stall cell was much larger than the impeller stall 
cell and strongly influenced the compressor performance 
during the off-design flow operation. 

(5) The evolution process of the LEV was precisely 
simulated by the in-house DES analysis and found to be 
the cause of the diffuser stall. A full annulus analysis 
pointed to the spillage of the LEV being related to the 
diffuser stall. 
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